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Power spectrum of mass and activity fluctuations in a sandpile 
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We consider a directed abelian sandpile on a strip of size 2 x n, driven by adding a grain randomly 
at the left boundary after every T time-steps. We establish the exact equivalence of the problem of 
mass fluctuations in the steady state and the number of zeroes in the ternary-base representation 
of the position of a random walker on a ring of size 3" . We find that while the fluctuations of mass 
have a power spectrum that varies as 1// for frequencies in the range 3~ 2n <C / <C 1/T, the activity 
fluctuations in the same frequency range have a power spectrum that is linear in /. 

PACS numbers: 



I. INTRODUCTION 

A subject of continuing interest in the study of non- 
equilibrium systems is the question of whether there is 
a common underlying mechanism in the very wide va- 
riety of systems where "1//" noise is observed [ll— fill] - 
More generally, the spectrum is a power law, 1// Q with 
a near 1. Further, the exponent a can be different for 
the noise in different physical quantities in the same sys- 
tem. In the past several decades it has been realized that 
1// Q spectra can arise in a number of different contexts, 
including fractal renewal processes fl2| . constrained re- 
versible Markov chains [13| , nonlinear stochastic differ- 
ential equations [14| . or multiplicative random processes 

Bak, Tang, and Wiesenfeld (BTW) proposed that for 
a large class of driven non-equilibrium systems, the 1/ f a 
noise is a manifestation of the long-ranged temporal cor- 
relations of fluctuations in their self-organized critical 
(SOC) steady state |1||[I3j an d introduced simple driven 
automaton model of sandpiles that reached a state char- 
acterized by power-law time and space correlations. Crit- 
icality does not necessarily imply that the exponent a of 
the power-spectrum is approximately 1 |18l420|. 

Studies of noise in self-organized critical systems have 
generally used the sandpile models as a paradigm. Most 
of these have been made in the slow driving limit, when 
the time-interval T between consecutive particle addi- 
tion events is much greater than the typical duration of 
an avalanche. In this limit, for frequencies 1 ^> / 3> 
1/T, the time-correlations in the noise-signal can be ex- 
pressed in terms of the critical exponents that charac- 
terize the scaling of the probability distribution of a sin- 
gle avalanche, and the exponent a is also expressible in 
terms of these exponents [21(. Numerical results indi- 
cate that the activity in the BTW model exhibits a non- 
trivial 1// Q power spectrum, with the exponent a be- 
ing between and 2 21 1. Some variants of the BTW 
model also exhibit a nontrivial 1// power spectrum [22| - 
[25j . Numerical simulations of the avalanche spectrum of 
the Manna model in 1, 2, and 3 dimensions gave expo- 



nents that ranged from 1.44 to 1.9 [21|. The 1// noise in 
the fluctuations of mass was first seen in a sandpile model 
with threshold dissipation (the sand is not conserved in 
topplings when the initial height is too large) by Ali [26| . 

For 1-dimensional models, scaling arguments and nu- 
merical studies have shown that in the single avalanche 
regime, the power-spectrum of the fluctuations in the to- 
tal activity has an exponent a=l. There is much less 
work dealing with noise correlations for durations much 
greater than T . These are governed by correlations be- 
tween different avalanches, which are much harder to de- 
termine. Hwa and Kardar [27| identified three different 
regimes of frequencies or time-scales, corresponding to 
typical duration of a single-avalanche, interval between 
of successive particle additions, and the relaxation time 
of the medium. They argued that the effective expo- 
nent a takes different values in different regimes. For 
/ <C 1/T , the avalanches are anticorrelated, and one typ- 
ically finds that the power spectrum varies as f a , with 
a' > 0. Maslov, Tang and Zhang (MTZ) HH studied a 
one-dimensional directed model of sandpiles, and showed 
that o;=l for noise in the spectrum of fluctuations of mass 
for / < 1/T. 

The above fact — and the analytical tractability of 
1-dimensional sandpile models — make the study of 1 
dimensional sandpile models especially interesting, and 
in the present paper we study the model of sandpiles in- 
troduced by MTZ in more detail. For a sandpile on a 2 x n 
strip, while the mass fluctuations have a 1// spectrum for 
a very wide range of frequencies 1 ;§> jT 3> 9 _n , the fluc- 
tuation of activity in the same range of frequencies show 
a power spectrum that is proportional to /. We show 
that mass fluctuations in this sandpile can be analyzed 
in terms of the position of an unbiased random walker 
on a 1-dimensional ring of 3™ sites. The mass of the pile 
is related to the number of zeroes in a ternary represen- 
tation of the position of the walker. While the random 
walk on a ring is a rather classical problem, the number 
of zeroes in the binary (or ternary) representation of its 
positional coordinate is an interesting object that does 
not seem to have been studied before. It is obvious that 



this has very slowly decaying correlations, and so long as 
the mean square deviation of the walker varies as a power 
of the time, the variance of the change in this number in 
time t is expected to vary as logt. This in turn implies 
that spectrum of mass fluctuations has a 1// tail for a 
very large range of frequency. Note that this mechanism 
underlying 1// spectrum differs from earlier approaches 
that obtain it as a superposition of multiple Lorentzians 
coming from a distribution of relaxation times [25| . 

The plan of this paper is as follows. In Section II we 
briefly describe the sandpile model on a 2 x n strip. We 
then present the details of the ternary base representa- 
tion, and the mapping to the random walk on a ring. In 
Section III wc calculate the time correlation of mass fluc- 
tuations in the sandpile using the new ternary represen- 
tation and in Section IV we compute the local activity, 
its mean, and the asymptotic velocity of an avalanche 
wave. In Section V we calculate the mean total activ- 
ity and present numerical results for the power spectrum 
of its fluctuations. The activity fluctuations scale as the 
time-derivative of mass and thus have a power spectrum 
proportional to /. Section VI contains a summary of our 
results. 



II. DEFINITION OF THE MODEL AND ITS 
RELATION TO A RANDOM WALK 




FIG. 1: The directed sandpile on a 2 x n strip (left) is equiv- 
alent to a random walk on a ring with N — 3™ sites (right). 
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The model studied earlier by MTZ is a directed sand- 
pile [29| on a strip of dimension 2 x n, shown in Fig. [T] 
An integer variable z(x,y) is assigned to each site (x,y) 
of the lattice: this is the "number of grains of sand" (or 
height) at site (x, y). The critical height is z c = 1 and if 
z(x,y) > z c , that site topples and transfers two grains, 
one to its right neighbor and one to the other site in the 
same column. One grain leaves the pile for each toppling 
that occurs at the rightmost column, x = n. 

The sandpile is driven periodically with a grain of sand 
added after every T time steps at one of the two sites in 
the first column, selected randomly. If as a result of this 
process the site becomes unstable, we relax the configu- 
ration by toppling all unstable sites in parallel, until all 
sites are stable. Further, if T > In the avalanche activity 
will have died out before the next particle is added [3fJ. 

As noted by MTZ [28[ the number of recurrent config- 
urations in this model is N = 3™. In a recurrent configu- 
ration, each column r can have the configuration ( J ), ( } ) 
or ( 1 ) , independent of the status of other columns. De- 
fine an integer Wj which takes values +1, or —1, if the 
configuration in the j th column is ( J ) , ( \ ) or ( \ ) respec- 
tively. Then to each of the 3 n recurrent configurations, 
we can attach a unique label j given by 



£>3 r 



(1) 



FIG. 2: (Colour online) Typical evolution of a sandpile with 
n — 4, seen as a random walk on a ring of size 81. 



counted from left to right as usual, and it is convenient 
to adopt the convention that the ternary representation 
of the integer j, namely r\\r\i ■ ■ ■ r\ n is also read left to 
right. Note that then the configuration with the label 
1010 corresponds to the integer j = 1 ■ 3° + ■ 3 1 - 1 • 3 2 + 
0-3 3 + c, withe =40 and not 0-3°- l-3 1 + 0-3 2 + l-3 3 + c. 
Let Aii and A\^. denote the operators correspond- 
ing to adding a particle at site (1,1) and (1,2) respec- 
tively, and relaxing. These satisfy the abelian property 
[Ai,i, A-1,2] = 0. Let the height configuration represented 
by the integer j be denoted by \j). Then wc find that 



Ai,i\j)->\j + iy, A x>2 \j)^\j-i) 



(2) 



where the addition is understood to be mod 3™, so that 
|3") = |0). 

The proof is quite straightforward. For any unstable 
configuration with heights {z(x,y)}, we define the vari- 
ableF({z( a; ,y)})by 



F({z(x,y)}) = y £ l 3<'- 1 [z{x,l)-z(x,2)]. 



(3) 



with c = (3" - l)/2, so that < j < 3" - 1. Rows are 



It is easy to sec that by adding a particle at (1,1) F 
increases by 1, while by adding at (1,2) it decreases by 



1. A toppling at any site (x, y), with i/b leaves F un- 
changed, and a toppling at the rightmost column changes 
F by ±3™. 

This can also be seen from the abelian algebra of the 
sandpile [31]], using the fact that Ai : i and A\ t i can be 
considered as generators of the Abelian group, with the 
relations 1281 



^1,1^1.2 = I 



N 



K 



N 



K 



(4) 



Representing the N recurrent configurations of the pile 
as sites on a ring (see Fig. [1} , stochastic addition of par- 
ticles to the pile at the left end gives a random walk: at 
each time step, the walker at site j has equal probability 
to take a step to the left or the right. We will denote the 
random walk corresponding to a particular realization of 
the random evolution of the pile by W, and the position 
of the walker at time t by X t . A typical realization of 
such a random walk is shown in Fig. [2] for a lattice of size 
n = 4. 



III. THE AUTOCORRELATION FUNCTION OF 
MASS FLUCTUATIONS 

The total number of grains in the pile at time t is the 
mass, denoted by m(t). The autocorrelation function 
C(t) of the fluctuations in mass in the steady state of 
the pile is 



C{t) = (m(t)m(t + r)) - (m(i)} 2 



(5) 



where the angular brackets denote averaging over the 
steady state. 

It is easy to see that (m(t)) = An/3. Let rrij(t) denote 
the mass of particles in the column j. Since in the steady 
state, the masses rrij in different columns are independent 
random variables, the equal time mass correlation is 



C(0)=£[(m|>-K> 2 ] 



(6) 



We define a variable £j corresponding to the jth bit of 
the position X t of the walker by £j(X t ) = 1 if r}j = 0, 



and £j(X t ) 
pile is 



otherwise. Then the mass m(t) of the 



m(t) 



3=1 



(7) 



The function £j (X) is a (non-random) periodic func- 
tion of its argument X of period Nj = 3 J . For example, 
£i(X) is 0,1 or — 1, for X= 0,1,2 (mod 3) respectively 
Expressing £j(X) in terms of its Fourier series expansion 



fc(J0=£6-(fc)exp(tfcr) 



(8) 
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FIG. 3: Part of the time-series of mass fluctuations £ for a 
system of size n = 6. 



where k = 2ttt/Nj and r = 0,1, ... ,Nj — 1, the coeffi- 
cients are given by 



N 



Nj-ik^ 



1 e-^-Olr g^HjzJ 
Nj [ sin (|) 



(9) 



Consider the auto-correlation function fj (r) of the ran- 
dom variables £j , defined as 

fj{T) - ^(X t )^(X t+T )) 

= ££&(fe)&(*0(e^ t+ife ' Xi+T >- (10) 

k k' 

By translational invariance along the ring, this expecta- 
tion value is zero unless k + k' = 0. Also, ( e lkx t-ikx t+T ^ 
is the characteristic function of the displacement X(t + 
r) — X(t) of a simple random walk W and this is easily 
calculated, giving 

/ 3 -(r)=El^(*)| 2 [cosfc]W. (11) 

k 

Using Eq. IJ5J), this can be evaluated, and the result is 

(12) 



f j (T)=Y,[ ^ 2 - ■ 1 2 [CC 



where k r = 2-irr/Nj, with r = 0,1, ... ,Nj — 1. For fc r < 
1, we can write [cos fc r ]' T w exp(— fc^|r|/2). Also, Nj-ik r 
only takes values 0, 2tt/3, 4tt/3. 

The r=0 term is the disconnected part of the corre- 
lation function, and this gives a finite contribution, 1/9. 
When r is a nonzero multiple of 3 the numerator van- 
ishes, and we have 



fj ( T )Kh(\T\/V) + l, 



(13) 
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FIG. 4: (Colour online) (a) Autocorrelation function Cj(t) 
of time series £,■ for different values of j for system size n — 
6. The discrete data points are joined by lines as an aid to 
the eye; the theoretical curves [Eq. (|12|l ] are indistinguishable 
from the numerical results, (b) Data collapse of Cj(r) for 
different j onto a single curve with scaling variable it = 9 _J r. 



where 



h ^ = Y!jZ^ e M~27T 2 r 2 u) 



(14) 



The prime denotes that r is not a multiple of 3. As a 
simple check one can verify that this gives h(0) = 2/9. 

Since £,■ and £,-' have different periods, these time series 
are uncorrelated, and in the steady state the average of 
the product $^,-_t 7 -/(£ i j(t)£j'(£ / )) reduces to a product of 
averages of the individual terms, and thus 






(15) 



Therefore the connected part of the correlation function 
can be reduced to 



C(t)=]TQ(t) 



i=i 



3=1 



h(9~ 



(16) 
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FIG. 5: (Colour online) (a) Plot of the autocorrelation func- 
tion of mass fluctuations for different system sizes n. (b) Col- 
lapse onto a single scaling function when the scaling variable 
u = 9~"r is used. 



In the range 9 s < r < 9 S+1 the first s terms in the above 
summation are nearly zero, and each of the remaining 
terms is nearly equal to h(0). Thus we have 



C(r)« (2/9) [n-tog 9 |r[ 



(17) 



Since the global correlation function varies as — logr, 
the resulting power spectrum S(f) will therefore show 
an asymptotic l/f a dependence with a = 1, for 1 3> 
/ > 9-". 

Results are shown in Fig. [3] for the time-series of mass 
fluctuations £. In Fig. 01 we have shown the observed 
correlation functions Cj(t) for different values of j. The 
data was obtained by a Monte Carlo simulation for sys- 
tem size n = 6, using a time-series of £,■ of total length 
10 8 and sample averaged over 10 2 realizations. These are 
also compared with the theoretical values of these func- 
tions for the corresponding correlation functions Cj(t) 
using Eq. (fl2"]l. 

In Fig. [5j we have also shown the auto-correlation 
function of fluctuations of the total mass for different 




FIG. 6: (Colour online) Power spectrum of mass fluctuations 
for different n showing the 1// dependence. 
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FIG. 8: (Colour online) The mean activity of a column at 
distance x for a system of size n = 10, and T = 25. The 
topmost curve (open circles) denotes the sum of the local 
activity in the different columns. 



values of n. We find that data for different n collapse onto 
a single scaling curve when the scaling variable u = t/9" 
is used. The power spectrum of these mass fluctuations 
is shown in Fig. [5] for different system sizes. 



IV. FLUCTUATIONS OF THE LOCAL 
ACTIVITY 



Define the local activity a(x, t) of site x = (x,y) at 
time t as 1 if the site x topples at time t and other- 
wise. Addition of a particle at the left edge results in an 
avalanche wave, the evolution of which can be graphically 
depicted by drawing a directed arrow along a bond if a 
particle was transferred along it during the avalanche. A 
typical such wave is shown in Fig. [7] The avalanche activ- 
ity can branch, and these branches may rejoin. The max- 
imum time that a wave can take to propagate through 
the lattice is clearly 2n. 

The average activity (a(x, t)) at the site x at time t in 
the steady state of the driven pile is a periodic function 
with period T. The mean local activity (a(x, t)) at site 
x = (1,1) can be computed easily. A particle is added at 
the top or bottom site with probability 1/2, and if added 
at (1, 1) and that site becomes active, then the site was 
originally occupied, the probability of which is 2/3. The 
mean local activity at (1,1) at time t = 1 is therefore 



5 • | = 1/3. if the particle is added at the bottom site 
and it becomes active, then at the next time step it will 
be active only if the left edge had the configuration ( \ ) , 
and this has probability 1/3. The avalanche wave begins 
with R B U (symbols R B =— >, Rt =~* , U =t and D =\, 
are used to show arrow representation of avalanche wave) 
and the mean local activity (a(x, £)) at site x = (1, 1) at 
t = 2 is 4 • i = 1/6. The top-bottom symmetry implies 
that the mean local activity of the bottom site is equal 
to the mean local activity of the top site. 

In a similar manner, the mean local activity (a(x, t)) 
at site x = (2,1) can be calculated. The minimum time 
taken for the avalanche to reach this site is 2. If the wave 
is RtRti it is simple to see that the mean local activity 
is i • (|) 2 = 2/9. The maximum time taken is 4 and 
the wave is RtDRbU giving activity \ ■ (|) 2 = 1/18. At 
time 3 both RbURt and RbRbU are responsible, giving 
activity 2 • \ ■ \ ■ § = 2/9. 

Generalizing this, one can easily deduce that the local 
activity of any site x = (a;, y) at time step t is (r = i 
mod x, i = t mod T and p = 1/3) 



(a(X)i)> = (1(^(1 -P)- **<i<2*; 
x v ' " 10 otherwise. 



(18) 



The time taken by the avalanche to reach site n is T n = 
n + r with n < T n < 2n. The asymptotic velocity of the 
avalanche can be computed, 



(v) = lim (— ; 

n— >oo J n 



3/4 



(19) 



FIG. 7: A typical avalanche wave with branches. 



since the mean local activity follows a binomial distribu- 
tion, and (r) — n/3. The mean local activity, shown in 
Fig. [5J has a binomial distribution. 
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FIG. 9: Time evolution of total activity A(i) in the pile as a 
function of time t for a system of size n = 7, and with T = 20. 
The activity A(t) takes values 0, 1, or 2. The mean activity 
A(t) and the excess activity 6A(t) are also shown. 



V. THE MEAN TOTAL ACTIVITY 

Let A(i) denote the total activity in the pile at time t, 
A(t) = Y,a(x,t). (20) 

X 

Its mean, averaged over all possible evolutions is A(t) and 
the fluctuations are SA(t) = A(i) — A(t). The power- 
spectrum of A(t) consists of two parts: 8- function peaks 
at integer multiples of frequency wo = 2tt/T coming from 
A{t) which is a periodic function of time with period 
T, and a continuous part coming from the fluctuations 
8A(t). 

We now discuss the exact calculation of A(t). This 
is important for numerical studies since the power- 
spectrum is dominated by the 8— function peaks, and only 
after subtracting the periodic part it is possible to obtain 
a good estimate for the continuous part of the noise- 
spectrum. 

Let Fj be the particle flow out of the jth column. 
There are five options: No flux (</>), T,B,TB 7 and BT, 
where T denotes a particle coming out on the top row, 
B out of bottom row, TB two particles, one from the 
bottom and then one out of the top row (note that the 
later event is on the left), and similarly, BT. 

The time evolution of F as the avalanche moves along 
the strip can be deduced in a straightforward manner. 
Given Fj and the configuration of the sandpilc at column 
\j + 1), Fj + i we can be determined using the toppling 
rules: 



T. 



T and the configuration is ( J ), then -Pj+i = 



BT. 



Define the ket vector \a) = \<j>),\T),\B),\BT), and 
\TB) for a = 0, 1, 2, 3, 4. Let Prob(j, a, r) be the proba- 
bility that the flux at column j is of type a, and starts 
after r time steps. The speed change in the avalanche 
wave can be accounted for by attaching a weight z at 
each time step, and defining 



|P,-) = ^z T Prob(i,a,r)|a). 



(21) 



The initial condition is 



\Po) = ^\T) + -\B) 



from which it can easily be worked out that 

|Pi) = ±\4>) + |[|T) + \B) + \TB) + \BT)]. 

Since each column can have only one of three configura- 
tions, (J ), ( $ ) or (i), the evolution of the toppling wave 
as it moves along the strip is governed by the equation 



\P j+1 ) = M\Pj) 



(22) 



where, M = |(Mi + M 2 + M 3 ) and M x is a matrix that 
specifies the change in | Pj ) when it encounters ( i ) . Ex- 
plicitly, 





/l 1 0\ 

z 


Mi = 







z z 2 
\0 0/ 


It can be easily verified that Mi\<fi) 


z\T),M x \TB) = z 2 \BT). Similarly, 


/l 1 0\ 




M 2 = 


2 





\0 z 2 z) 



= |0),Mi|T) = 



and 



Mo = 



/l 0\ 

ooooo 1 


z z 

\0 z z/ 



If Fj = F and the configuration is (J), thenFj + i = Then, the mean activity A(t) at time t is given by the 



coefficient of z* in X^i^l-Pj) where 



(0,1,1,1 + 




o 





T=16 

T = 24 

T = 32 


1 1 


1 1 


4 -3 
Logf 


-2 -1 



matrix. \Pj) is always of the form a 3 -(\T) + \B))+/3j \TB) 




FIG. 10: (Colour online) The power spectrum of the fluc- 
tuations of excess activity 8A(t) in a sandpile of size n—7. 
(a) The spectrum for different T using a time series of total 
length 5 x 10 6 showing the initial linear portion as well as the 
1// part, and (b) the spectrum for the case when T is not 
fixed (a new particle is added when there is no activity)from 
a time series of total length 10 6 . The power spectrum has 
been binned using 5 bins per decade of frequency, and a line 
of slope -1 is shown for comparison. 



z, 1 + z\ and the Fourier or Laplace transform of A(i) is 



A(z) = £>|i>-> 



±m Ml+ T M * *™- ^ 



A(z) = J2( 1 l + «) 





z(z+2) 



A(z) 



EC 11 



3 

Z 22 

3 3 



(24) 



j'=i 



and this tridiagonal matrix is easy to evaluate. We get 
A(z) = 2z/3 + 7z 2 /9 + . . . . The mean total activity at 
time step k can be written as 

1 - (1/3) + (1/3 2 ) • • • + (1/3*) = |[1 - (-i) fc+1 ] (25) 

which goes to 3/4 as k — ► oo. 

Shown in Fig. [9] is a typical realization of the total ac- 
tivity A(t) for n = 7. The corresponding power spectrum 
is shown in Fig. [TUJ We see that there is a weak depen- 
dence of the spectrum on T, the period of addition of 
particles, shown in Fig. HPT a) and the power spectrum of 
these signals exhibit asymptotic / behavior for / <C 1/T. 
This behavior may be understood as follows: the fluctu- 
ations of activity scale the same way as the fluctuations 
of number of particles leaving the system. When this 
number is non-zero, the corresponding avalanche spans 
the system, and the total activity is also large. When 
the avalanche wave does not span the system, the activ- 
ity is less, and the outflux of particles is zero. Hence, 
we expect that scaling properties of time-correlations of 
fluctuations of total activity to be similar to that of fluc- 
tuations of outflux. But the time series of the latter is 
obtained by taking successive differences of time-series 
for fluctuations of mass. Then, the power spectrum for 
the fluctuations of total activity is related to that of mass 
by an additional factor f 2 . 

Interestingly, if the time interval between additions of 
grains is not fixed and a new grain is added as soon 
as the avalanche generated by the previous grain has 
stopped, the behavior of activity correlations changes. 
In Fig. II Or b), we have shown the result for power spec- 
trum of total activity fluctuations in such a simulation. 
We find that the power spectrum shows 1// behavior for 
several decades of frequency range in this case. It has 
been noted earlier that the power spectrum for frequen- 
cies / > 1/T is not substantially affected by this change 
[21j . but the reason why this extends the range of 1// 
scaling to much lower frequencies is not understood. 



VI. SUMMARY AND DISCUSSION 



The matrix calculation can be simplified by first noting 
that \4>) can be dropped from the calculation since it does 
not contribute. Then from the symmetry of the top and 
bottom rows the 4x4 matrix can be reduced to a 2 x 2 



In this paper we have established an exact equivalence 
between the total mass-fluctuations in a directed sand- 
pile on a ladder and the number of zeroes in the ternary 
representation of a random walker on a ring. We obtained 
an exact expression for the time correlation function of 



s 



mass fluctuations. Local temporal correlations have a 
characteristic exponential dependence on the spatial co- 
ordinate and from an exact calculation we find that the 
functional form of the global correlation is logarithmic, 
resulting in an 1/f spectrum for frequencies below a cut- 
off i / < fc- We also calculated the mean local activity, 
the mean total activity and the asymptotic velocity of 
this avalanche wave was shown to be (v) = |. For the 
periodically driven case the power spectrum of total ac- 
tivity fluctuations exhibits 1/f dependence for frequen- 
cies / satisfying 1 > /> 1/T, and a power proportional 
to / for 1/T > / > 3~ 2 ". When a new particle is added 
as soon as an avalanche is finished the power spectrum 
is 1/f over the entire range l>/> 3~ 2n . 

It is also important to note that in this model the lower 
cutoff on the range of observed power-law dependence is 
very small, varying as exp(— L) where L is the linear ex- 
tent of the system. There are related models wherein 
events relax very very slowly: for instance in ricepiles 
where the toppling condition depends on the slope, the 
mean residence time for a marked grain to leave the pile 



varies as cxp(aL d+2 ) for a d-dimensional system [32JJ. 
It would be interesting therefore to investigate how the 
lower cutoff for the 1/ f a spectrum depends on system 
size in other sandpile models, particularly when there 
are properties that can have exceedingly long relaxation 
times. 
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